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Abstract 

Ideally,  an  adaptive  optical  control  system  would  have  instantaneous  temporal  response  and  infinite  spatial 
bandwidth.  In  real  systems,  the  response  time  of  the  adaptive  optical  control  system  is  limited  by  the  integration 
time  of  the  wave  front  sensor,  the  computational  time  of  the  control  algorithm,  and  the  actuator  response  time. 
Additionally,  finite  inter  actuator  spacing  limits  the  deformable  mirror’s  ability  to  reproduce  spatial  frequencies 
having  a  period  less  than  twice  this  spacing.  Although  analyses  general  enough  to  account  for  both  the  temporal 
and  spatial  characteristics  of  the  adaptive  optical  system  exist,  they  are  complex  and  require  detailed  information 
regarding  the  wave  front  sensor,  the  deformable  mirror,  and  the  control  algorithm.  Ibis  investigation  develops 
a  frequency  domain  model  that  describes  performance  effects  of  an  adaptive  optical  system’s  temporal  response 
taking  into  account  aperture  piston  and  tilt  removal  and  spatial  bandwidth  limitations  due  to  finite  subaperture 
size.  The  unique  aspect  of  this  model  is  the  relative  ease  with  which  performance  characteristics  of  different 
spatial  and  temporal  system  response  functions  can  be  investigated. 


The  Effect  of  an  Adaptive  Optical  System’s  Spatio-Temporal 
Response  on  Imaging  Performance 


I.  Introduction. 


1.1  Background. 

Adaptive  optics  are  used  to  improve  the  real  time  imaging  capabilities  of  large  ground  based  telescopes 
(9).  As  light  from  a  celestial  source  passes  through  the  earth's  atmosphere,  it  suffers  random  phase  delays  across 
its  wave  front.  These  phase  delays  are  caused  by  turbulence  induced  fluctuations  in  the  atmosphere’s  index 
of  refraction.  Because  of  the  effect  of  atmospheric  turbulence,  very  large  diameter  telescopes  have  little  more 
resolving  power  than  telescopes  0.1  to  0.2  meters  in  diameter  (15)(19).  By  measuring  and  compensating  for 
wave  front  distortions,  adaptive  optics  make  it  possible  to  achieve  near  diffraction  limited  performance  of  these 
same  very  large  telescopes  (4).  Figure  1.1  depicts  a  typical  adaptive  optic  system.  First,  turbulence  induced 
phase  distortions  of  a  reference  wave  front  are  measured  by  a  wave  front  sensor.  This  phase  information  is 
then  used  by  a  computer  to  determine  the  optimal  actuator  weightings  across  the  deformable  mirror.  Typically, 
the  computer  implements  an  algorithm  that  minimir.es  the  mean  square  error  between  the  unabenrated  reference 
wave  front  and  the  wave  front  predicted  to  anive  at  the  image  sensor.  Lasdy,  the  phase  distortions  are  corrected 
at  the  deformable  mirror  surface  by  complementary  path  length  adjustments  across  the  wave  front,  and  the  object 
of  interest  is  imaged. 

Ideally,  an  adaptive  optical  control  system  would  have  instantaneous  temporal  response  and  infinite  spatial 
bandwidth.  In  real  systems,  the  response  time  of  the  adaptive  optical  control  system  is  limited  by  the  integration 
time  of  the  wave  front  sensor,  the  computational  time  of  the  control  algorithm,  and  the  actuator  response  time. 
Additionally,  finite  inter  actuator  spacing  limits  the  deformable  minor's  ability  to  reproduce  spatial  frequencies 
having  a  period  less  than  twice  this  spacing  (19).  As  any  one  of  these  system  parameters  can  be  the  limiting 
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Figure  1.1.  Representation  of  an  adaptive  optical  system  correcting  for  atmospheric  tubulence. 


factor  in  the  system’s  imaging  performance,  it  is  invaluable  in  system  design  to  have  simple  equations  or  models 
that  provide  a  point  of  departure  for  each  parameter. 

Response  time,  from  wave  front  detection  through  actuator  positioning,  limits  system  performance  because 
the  phase  distortions  across  the  wave  front  are  continually  changing  in  time.  This  results  in  an  incorrect  matching 
between  the  deformable  minor  and  the  aberrated  wave  front,  especially  for  higher  frequencies.  In  two  recent 
papers.  Fried  and  Kan  independently  investigated  the  effect  of  time  delay  on  the  servo  control  loop  of  an  adaptive 
optical  system  (11)(2).  Although  their  approaches  are  slightly  different,  both  made  the  same  assumptions  and 
determined  the  mean  square  phase  error  between  the  actual  wave  front  and  the  sensed  wave  front  to  be: 

<t2  =  28.4(/cAt)S  (1.1) 


where  fa  is  the  well  known  Greenwood  frequency  and  At  is  time  delay  between  wave  front  sensing  and  actuator 
repositioning  (6).  In  their  derivations  Fried  and  Karr  neglected  the  spatial  frequency  limitations  of  the  deformable 
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mirror  and  approximated  the  atmospheric  turbulence  power  spectrum  by  its  high  frequency  asymptote  ignoring 
the  inherent  removal  of  piston  and  the  possibility  of  tilt  correction  over  the  aperature.  They  also  neglected 
temporal  bandwidth  limitations  due  to  the  finite  integration  time  of  the  wave  front  sensor,  causing  an  unrealistic 
equal  weighting  of  all  temporal  frequencies  across  the  spectrum.  By  equally  weighting  all  frequencies.  Fried  and 
Karr’s  results  also  reflect  a  higher  than  expected  value  of  the  mean  square  error  for  any  time  delay  At.  As  the 
delay  time  becomes  longer,  lower  spatial  frequency  corrections  of  the  deformable  mirror  become  uncorrelated 
with  the  incoming  wave  front  increasing  the  residual  mean  square  error. 

The  above  cited  analyses  have  broad  applicability  to  systems  that  have  both  a  wide  spatial  bandwidth 
and  very  fast  response  time.  Adaptive  optical  systems  that  meet  these  requirements  are  normally  built  only 
for  smaller  aperture  telescopes.  The  complexity  and  expense  involved  in  producing  adaptive  optical  systems 
that  attempt  to  meet  these  requirements  for  larger  telescopes  have  driven  designers  to  considering  partial  spatial 
compensation  (16).  Systems  that  implement  partial  compensation  have  a  lower  spatial  bandwidth  and  are  not 
adequately  modeled  by  Eq  (1.1).  Although  analyses  general  enough  to  account  for  both  the  temporal  and  spatial 
characteristics  of  the  adaptive  optical  system  exist,  they  are  complex  and  require  detailed  information  regarding 
the  wave  front  sensor,  the  deformable  mirror,  and  the  control  algorithm  (22)(21). 

1.2  Objective 

The  primary  objective  of  this  investigation  is  to  develop  a  frequency  domain  approach,  using  linear 
systems  methods,  that  describes  performance  effects  of  an  adaptive  optical  system's  temporal  response  taking 
into  account  aperture  piston  and  tilt  removal  and  spatial  bandwidth  limitations  due  to  finite  subaperture  size.  The 
secondary  objective  is  to  address  when  use  of  our  approach  yields  substantially  different  results  than  those  of 
Fried  and  Karr. 
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1.3  Approach 

In  this  investigation,  a  frequency  domain  approach  is  used  to  predict  the  mean  square  value  of  the  residual 
phase  error  of  a  corrected  wave  front,  ipc(x,t).  Modeling  the  system  transfer  function  of  the  adaptive  optical 
system  as  shown  in  Figure  1 .2,  the  residual  phase  error  of  tpc(x,t)  can  be  shown  to  equal: 


<r2=  f°°  df  \\-H{f)?W{f) 

J -CO 


(1.2) 


where  W(f)  is  the  temporal  power  spectrum  of  the  wavefront  phase  fluctuations  due  to  atmospheric  turbulence. 
p(x,t),  and  //(/)  is  the  transer  function  of  the  adaptive  optical  system.  The  input  spectrum,  W(f),  is  derived 
as  the  Fourier  transform  of  the  temporal  covariance,  C^(t)  =  (<p(x,t)<p(x,  t  +  r)},  and  will  be  developed 
to  account  for  inherent  piston  removal  and  tilt  removal  capability  of  present  adaptive  optical  systems  using 
methods  developed  by  Greenwood  and  Fried  (8).  The  filtering  action  of  the  transfer  function  will  be  assumed 
to  be  separable  in  time  and  space  so  that  H(f)  =  Hj{f)  ■  Hsif  ).  where  HT{f)  accounts  for  the  system’s 
temporal  delay  and  bandwidth  and  Hs(f)  is  the  temporal  equivalent  spatial  response  of  the  deformable  mirror. 
It  is  important  to  note  here  that  all  filtering,  both  temporal  and  spatial,  will  be  done  in  the  temporal  frequency 
domain.  For  a  given  wave  front  spectra  m,  W  ( / ) ,  this  approach  lends  itself  to  the  investigation  of  arbitrary  filter 
functions  due  to  the  simple  multiplication  operations  required  to  account  for  the  filtering  response  in  Eq  (1.2). 
This  will  be  shown  to  be  especially  true  of  the  temporal  transfer  function 
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1.4  Overview  of  Thesis 


Chapter  II  reviews  Greenwood  and  Fried's  method  of  developing  a  piston  and  tilt  removed  spectrum  of 
wave  front  corrections  for  an  adaptive  optical  system.  Chapter  HI  developes  the  temporal  spectrum  of  wave  front 
phase  aberations  a  counting  for  piston  and  tilt  and  developes  a  method  for  determining  the  temporal  equivalent 
spatial  response  of  the  deformable  mirror.  In  chapter  IV  the  result  of  the  preceding  chapter  is  evaluated  for  values 
of  interest  to  current  system  designers  and  compared  to  results  using  Karr’s  and  Fried’s  models.  Conclusions  are 
drawn  and  recommendations  made  in  chapter  V. 
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II.  Literature  Review. 


This  chapter  reviews  Greenwood  and  Fried’s  (GF)  often  referenced  paper  which  describes  the  power 
spectra  of  wave  front  corrector  motion  for  circular  correcting  subapertures  of  a  larger  segmented  circular  mirror 
(8).  In  their  analyses,  the  correcting  subapertures  of  the  system  are  capable  of  independent  piston  and  tilt  motion, 
and  the  power  spectra  for  both  types  of  motion  are  developed.  Additionly.  GF  account  for  overall  piston  and 
the  possibility  of  overall  tilt  correction  prior  to  the  wave  front  reaching  the  segmented  mirror.  It  is  this  aspect 
of  their  development  that  is  particularly  useful  in  this  investigation.  Figure  2.1  depicts  the  wave  front  correction 
process  modeled  by  GF.  Tilt  of  the  wave  front  is  removed  as  the  best  fit  plane  through  the  center  of  an  aperture. 
The  center  position  of  a  correcting  subaperture  is  given  by  r\  and  is  referenced  from  the  center  of  the  overall 
aperture,  (see  Figure  2.2).  By  allowing  the  subaperture  size  of  the  segmented  minor  to  approach  zero,  the  piston 
and  tilt  removed  spectrum  of  the  incident  wave  front  phase  can  be  approximated  by  the  piston  and  tilt  removed 
spectrum  of  the  subaperture’s  piston  motion  at  each  point  on  the  mirror. 

To  begin  statistical  modeling  of  a  subapertures  piston  motion,  it  is  first  necessary  to  have  an  expression 
for  the  vertical  position  of  the  subapertures  as  a  function  of  time.  GF  find  the  time  dependant  piston  position  of 
a  subaperture  of  diameter  d  centered  at  position  rj (  to  be 

<p(t)  =  -  <*?(£,*)  -  1 6nD~2fi  xip(x,t)  (2.1) 

where  a  single  overbar  denotes  a  spatial  average  over  the  subaperture,  the  double  overbar  denotes  a  spatial 
average  over  the  overall  aperture  of  radius  D,  and  ^p(x.t)  the  phase  of  the  wave  front  prior  to  any  corrections. 
The  first  term  of  Eq  (2.1)  represents  the  average  phase  of  the  wave  front  over  a  subaperture  centered  at  r\  in 
the  aperture  at  time  t .  Term  two  represents  the  average  phase  of  the  wave  front  over  the  entire  aperture  and  its 
inclusion  accounts  for  the  system's  insensitivity  to  aperture  piston.  Lastly,  term  three  accounts  for  the  possibility 
of  aperture  tilt  removal  by  the  adaptive  optical  system.  In  this  equation,  n  =  1  indicates  aperture  tilt  removal 
and  n  =  0  that  no  such  compensation  has  occured. 
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(  C  ) 


Figure  2. 1 . 


"A  one-dimensional  representation  of  wave-front  correction  (vertical  scale  greatly  expanded).  (A) 
Initial  wave  front  at  aperture  plane,  (B)  wave  front  after  overall  tilt  correction;  local  piston  and  tilt 
shown  for  segment  of  size  d  =  D/7,  (C)  residual  wave  front  after  compensation  over  all  segments. 
(8)." 
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aperture 
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Figure  2.2.  Top  view  of  a  subaperture  centered  at  position  f\ . 


To  determine  the  power  spectral  density  or  spectrum  of  a  real,  stationary,  zero  mean,  random  process,  we 
first  need  to  determine  its  covariance  function.  Fourier  transforming  the  covariance  function  results  in  the  power 
spectral  density.  The  temporal  covariance  of  <p(t )  is  defined  as 

(t)  =  (<p(t)<p(t  +  r)).  (2.2) 

GF  substitute  Eq  (2.1)  in  Eq  (2.2)  and  find  the  temporal  covariance  can  be  written  in  the  following  form: 

C*  (r)  =  J  dr  Ctip  (r,  r)  Tv  (r)  (2.3) 

where  (r,  r)  =  (ip  ( x ,  t)  ip  (x  +  f,  t  +  r))  and  Tv{r)  is  a  transfer  function  defining  the  effects  of  piston 
and  tilt  removal  at  single  points  in  the  aperture.  GF  find  (r)  to  have  the  following  form  for  a  subaperture  size 
that  approaches  zero: 

T,  (r)  =  6(r)  +  5j  (r)  -  S2(r)  (2.4) 


where 


Si  (rl  = 


{4«  (55-)  [(5J-)  +  (*)  cos(0 -<*)]  +  !} 

X  jarccos(£)  -  (£)  [l  -  (£)2]  'J 

-  |n  (if)  (%)2  [l  -  (^)2]  f  [1  +  4  cos’(0  -  *]  ,  0  <  (fr)  <  1 

0.  1  <  (b) 


(2.5) 


and 


7 tD2 


1+4"  (&)[(&) +  (£)«»(*-*].  0<(£)<i 

o.  *<(£) 


(2.6) 


and  r  =  rj  is  the  vector  specifying  the  point  of  interest  in  the  aperture,  6  =  arg  (r).  4>  =  arg  (fi ),  and 
r2  =  \Jrl  -f  r,  +  2  rri  cos  (6  -  <t>).  The  actual  derivation  of  the  transfer  function,  7^  (r),  ispresenredbyGF.lt 
is  important  to  recognize  at  this  point  that  (r )  represents  the  covariance  over  a  subaperture  with  zero  diameter 
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located  at  ?i  and  not  over  the  entire  aperture.  Recall,  we  are  considering  a  subaperture  size  that  approaches  zero 
diameter  to  approximate  the  incident  wave  front  phase. 

GF  use  the  following  relationship  to  determine  C'tv  (r ,  r): 

Ctlf,  (?,  r)  =  C6ifi  (0,  r)  -  (?,  r)  (2.7) 

where  (f,  r)  is  the  called  the  phase  structure  function  by  GF  and  is  defined  as: 

Df,^  ( r ;  r)  =  ([<p(x,t)  -<f(x+  r,t)][<p(x,t  +  r)  -  <p(x  -f  r,t  +  r)])  (2.8) 


Substituting  Eq  (2.7)  in  Eq  (2.3),  GF  anive  at  the  important  result: 

ct(T)  =  -1-Jdr  D6v 


(2.9) 


where  as  shown  by  GF: 

J  dr  Tv  (?)  =  0. 


(2.10) 


Note  the  inclusion  of  the  delta  function  term  in  Eq  (2.4).  GF  validate  Eq  (2. 10)  for  a  nonzero  subaperture  diameter. 
d.  and  when  discussing  the  single  point  phase  spectrum,  state  the  delta  function  term  can  be  disregarded.  This  is 
only  half  true.  The  delta  function  has  no  impact  on  Eq  (2.9)  but  is  in  fact  necessary  for  Eq  (2.10)  to  be  valid. 

Eq  (2.9)  can  be  Fourier  transformed  to  give  the  single  point  temporal  power  spectrum  of  piston  motion: 


=  dr  (r,  f)  (?) 


(2.11) 


where  IVif  (?,  /)  is  the  temporal  Fourier  transform  of  the  phase  difference  structure  function  introduced  in  Eq 
(2.8).  This  spectra  will  subsequently  be  referred  to  as  the  phase-difference  spectrum. 
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To  evaluate  Eq  (2.11),  only  (r,  /)  requires  definition.  GF  mistakenly  reference  Lee  and  Harp  (12) 
(they  should  have  referenced  Tatarski  (18) )  and  define  the  single  layer  phase-difference  power  spectrum  for  a 
plane  wave  reference  as: 

W'^(r,/)  =  2.079 (///n)-8^/*)-1  x  (D/r0)5/3  sin2(/r//„ )  (2.12) 

where  /„  =  57  is  a  normalizing  variable  that  accounts  for  wind  velocity,  v,  and  aperture  diameter.  D.  rD 
is  Fried’s  coherence  length  (a  measure  of  the  strength  of  a  layer  of  atmospheric  turbulence)  (3),  and  r  =  |f|. 
GF  arrive  at  Eq  (2.12)  by  simplifying  the  vectoral  nature  of  ( f,  /).  GF  assume  the  vector  r  and  the  wind 
velocity  v  are  parallel  for  every  pair  of  points  across  the  aperture.  The  significance  of  this  assumption  will  be 
clarified  in  the  next  chapter.  This  definition  of  the  phase  difference  power  spectrum  also  assumes  a  Kolmogorov 
spectrum  of  the  atmospheric  refractive  index  fluctuations  with  an  infinite  outer  scale  and  an  infinitely  small  inner 
scale. 

GF  evaluate  Eq  (2.1 1)  for  the  above  phase  difference  spectrum  and  present  their  result  in  graphical  form. 
For  perfect  spatial  wave  front  reconstruction,  d  =  0.  they  find  the  high  frequencies  roll  off  as  /-8/3  and  are 
unaffected  by  piston  and  tilt  correction.  Regardless  of  subaperture  size,  GF  find  the  low  frequencies  roll  off 
as  either  /-4/3  for  overall  tilt  removal  or  f~2/3  for  no  tilt  removal.  GF  also  provide  the  necessary  transfer 
functions  to  evaluate  Eq  (2.1 1 )  for  the  case  of  non-zero  subaperture  size.  The  resulting  spectra  would  include  the 
effects  of  piston  and  tilt  removal  and  also  finite  bandwidth  limitations  of  the  deformable  mirror.  Initially,  it  seems 
it  may  be  possible  to  use  the  asymtotic  approximations  to  represent  the  temporal  phase  spectrum  of  an  incoming 
aberrated  wave  front  corrected  for  overall  piston  and  tilt.  However,  the  GF  spectrum  plots  and  approximations 
are  not  valid  models  for  the  low  frequency  behavior  of  a  piston  and  tilt  phase  spectrum.  This  is  attributable  to 
the  non-realizable  assumption  that  the  phase  difference  power  spectrum  is  non-vectoral  in  nature  (7).  Eq  (2.1 1) 
and  the  transfer  function  ( r)  are  valid  however,  and  in  Chapter  HI,  these  equations  will  be  used  to  determine 
a  more  realistic  single  point  phase  power  spectrum  that  accounts  for  aperture  piston  and  tilt  removal. 
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111.  Methodology 


As  stated  briefly  in  Chapter  I,  Eq  (1.2)  will  be  used  to  model  the  mean  square  phase  error  performance 
of  an  adaptive  optical  system.  Before  this  equation  can  be  used  however,  an  expression  must  be  developed  for 
each  of  the  individual  terms  with  respect  to  the  temporal  frequency  variable  /.  Section  3.1  will  introduce  the 
spectrum  of  the  refractive  index  fluctuations  in  the  atmosphere  and  develop  the  associated  temporal  spectrum. 
Wv(f).  of  the  phase  fluctuations  <p(  x ,  t ) .  Section  3 .2  will  develop  the  aperture  piston  and  tilt  removed  temporal 
spectrum,  W$(f).  of  these  same  phase  fluctuations.  Lastly,  section  3.3  will  develop  the  method  by  which  spatial 
bandwidth  limitations  will  be  accounted  for  in  the  transfer  function  Hs(f). 


3.1  Phase  fluctuations  for  the  von  Karmen  spectrum 

This  section  introduces  the  von  Karman  spectrum  to  characterize  the  refractive  index  fluctuations  of 
turbulent  eddies  in  the  atmosphere.  From  this  spectrum,  the  spatial  and  temporal  spectra  of  the  wave  front  phase 
is  determined.  The  temporal  phase  spectrum  will  be  that  of  the  reference  wave  front  just  prior  to  aperture  tilt 
correction. 

The  3-dimensional  von  Karmen  spectrum  of  refractive  index  fluctuation  has  the  following  form: 


(k,)  = 


(3.1) 


where  k,  is  the  spatial  wave  number,  k0  =  ,  lcm  =  77,  and  C%  is  a  measure  of  the  strengh  of  the  atmospheric 

tubulence.  L0  and  l0  are  the  outer  and  inner  scales  of  the  atmospheric  turbulence  eddy  size.  This  spectrum  is  a 
variation  of  the  Kolmogorov  spectrum,  given  by: 


*n  (*t)  =  0.033C**,’"  .  (3.2) 

Goodman  notes  that  the  power  in  the  spectrum  below  k0  has  not  been  characterized  to  date,  and  that  the  von 
Karmen  spectrum  is  only  an  artificial  means  of  eliminating  the  pole  at  kt  =  0  in  the  Kolmogorov  spectrum  (5). 
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Given  $n  ( kt ).  we  can  detennine  Che  2-dimensional  spatial  spectrum  of  the  phase  fluctuations  of  a  plane 


wave  passing  through  a  turbulent  medium  as 


F  (kt)  =  nk2L 


*»(*«) 


(3.3) 


where  k  is  the  wave  number  of  the  incident  radiation,  and  L  is  the  length  of  the  propagation  path  through  the 
turbulence  (5).  Taking  advantage  of  the  rapid  decay  of  the  von  Kantian  spectrum  for  kt  >  km.  Eq  (3.3)  can  be 
approximated  as 

(kt)  w  2wPL$n  (kt)  (3.4) 


where  (kt)  is  defined  as  in  Eq  (3.1)  and  ¥«>■  This  assumption  is  satisfied  for  wavelengths  near  the 
visible  spectrum  and  is  typical  of  papers  in  this  research  area  (5)(18)(8).  Note  that  F(kt )  defined  in  Eq  (3.3)  is 
different  from  F^(kt)  defined  in  Eq  (3.4).  We  use  the  v  subscript  to  indicate  a  phase  spectrum  derived  from  Eq 
(3.1),  the  von  Karmen  spectrum. 

To  develop  an  equivalent  temporal  spectrum,  we  assume  that  the  structure  of  the  inhomogeneities  in  the 
atmosphere  remains  fixed  as  it  moves  with  the  local  wind  velocity  v.  The  wind  velocity  is  assumed  to  be 
perpendicular  to  the  propagation  direction  of  the  wave  front  (Taylor’s  hypothesis  of  frozen  turbulence).  The 
single-sided  temporal  power  spectrum  of  the  incident  wave  front  phase,  VV(/),  can  be  determined  from: 


W(f)  =  - 
v 


(3.5) 


where  the  derivation  of  Eq  (3.5)  is  presented  in  Appendix  A.  Substituting  Eq  (3.4)  into  the  above  expression: 


”V(/)  =  V 


mZCllnPL 


(3.6) 
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Letting  kf  =  +  ^o)  x  Eq  (3.6)  can  be  rewritten  as 


W'„  (/)  =  • 033  (Ji'j  *  C’Piv*  + 


p  n2 


(3.7) 


where  D  —  .  To  make,  this  function  integrable  in  closed  form  it  is  necessary  to  omit  the  exponential  term. 

This  can  be  done  if  is  assumed  to  be  large  enough  that  the  added  high  frequency  power  is  small  due  to  the 
remaining  polynomial  decay.  Using  the  Beta  function  identity 


(x  +  l)a+i 


T(a)T(b) 

T  (a  +  6) 


the  single-sided  temporal  power  spectrum  is  finally  written  as 


Wv(f)  =  .033C*k2Lv* 


(3.8) 


(3.9) 


3.2  Piston  and  tilt  removed  temporal  power  spectrum 

In  the  previous  section,  we  defined  the  temporal  spectrum  of  phase  fluctuations  for  a  distorted  wave  front. 
Here,  we  go  one  step  further  and  remove  the  power  in  the  wave  front  associated  with  aperture  piston  and  tilt. 
This  is  necessary  because  adaptive  optical  systems  are  insensitive  to  aperture  piston.  Additionally,  most  adaptive 
optical  systems  are  designed  to  remove  aperture  tilt  which  accounts  for  between  80%  and  90%  of  the  variance 
of  the  distorted  wave  front  (19X14).  By  removing  aperture  piston  and  tilt  from  the  input  spectrum,  we  can 
better  predict  how  well  an  adaptive  optical  system  will  correct  higher  order  aberations  as  well  as  more  accurately 
predict  temporal  response  requirements. 

The  removal  of  aperture  piston  and  tilt  can  be  accomplished  by  using  Eq  (2.1 1 )  and  the  transfer  function 
T,p  ( r)  developed  by  Greenwood  and  Fried  (GF).  In  Oder  touseEq(2.11)  and  avoid  the  problems  associated  with 
the  use  of  the  GF  approximations  to  the  phase  difference  power  spectrum,  W'^r,  /),  discussed  previously,  it 
is  necessary  to  first  develop  an  expression  for  W^(r,  f)  that  accounts  for  the  vectoral  nature  of  r.  The  vectoral 
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nature  of  Wiv( r,  f)  can  be  accounted  for  by  using  an  approach  very  similiar  to  that  used  by  Tatarski  to  derive  a 
relation  equivalent  to  Eq  (3.5)  (18). 

Following  Tatarski ’s  derivation, we  start  by  writing  the  phase  difference  structure  function: 

(r’»0  =  <[*>(*,<)  -  +  r,<)}  x  [<p(x-  vr,t)  —  <p(x  +  r  —  nr,  <)])  (3.10) 

where  the  two  terms  in  brackets  are  the  phase  differences  at  two  points  separated  by  distance  r  at  times  t  and 
t  +  r  and  a  frozen  flow  assumption  has  been  made  (18).  Tatarski  goes  on  to  show  that  the  identity 

(a  -  b)  (c  -  d)  =  1  [(a  -  df  +  (6  -  c)2  -  (a  -  cf  -  (6  -  d)2]  (3.11) 

can  be  used  to  convert  Eq  (3.10)  to  the  more  useful  form: 

Dio  (r,  r)  =  ^  {Do  (r  -  vr)  +  Do  (r  +  vr)  -  2 Do  (vt)}  .  (3.12) 

where  for  a  locally  isotropic  field  Dj  (r)  =  ([/  (r  +  r*)  —  /  (r’)]2). 

Utilizing  the  Fourier  transform  relationship  between  Dio  (r-  0  and  Wio  (/,  r). 


df  cos  ( 2 7r f  r )  Wio  (/•  r)  > 


and  the  relationship  between  Do  (r,  nr),  and  Fo(ht).  given  by  Lee  and  Harp  (12). 


Do  (r,  vt)  =  4tt  f  dkt  kt  [1  -  J0  ( kt  |r  -  nr|)]  Fo  (kt) 
Jo 


(3.13) 


(3.14) 
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and  finally  Eq  (3.12),  a  relationship  between  W6v>(r,  f)  and  Fv  (kt)  can  be  developed.  The  relation  can  be 
expressed  as  follows: 


2  Jo  cos(2,r^,r) =  JQ  d^< k*  t1  -  j°  (*«  \r-  ^rD]  (*i) 

+  4rr  [  d kt  k,  [1  -  J0  (kt  \f+  tTr|)]  F#  (kt) 

Jo 

-  8?r  J  dkt  kt  [1  -  J0  (£<t>r)]  F#  (fc,)j  (3.15) 


where  D^(x)  =  (x)  has  been  assumed.  Assigning  the  wind  velocity  a  particular  direction,  say  the  positive 

x  direction,  results  in  |r  —  vt\  —  \/r 2  +  v2t2  —  2rurcos$  and \r  +  tTr|  =  \/r2  +  v2r2  +  2rvr  cos  9,  where 
0  is  the  angle  between  rand  the  x-axis.  Eq  (3.15)  can  be  simplified  and  expressed  in  the  following  form: 

r  oo  roo 

2/  d/'  cos(2jr/'r)  Witfi(f',r)  =  2 ir  dktktF^(k,)x 

Jo  Jo 

[2 Jo  ( ktvT )  -  Jo  ( k ,  |r  —  nr|)  -  Ja  (kt  |f  +  nr|)]  (3.16) 

Eq  (3. 16)  can  be  further  simplified  by  Fourier  transforming  it  with  respect  to  t .  Using  the  following  relationship: 

/°°  1 

dr  cos  (2*/'r)  =  -[6  (/  +  /')  +  5  (/-/')]  (3.17) 

■OO  ^ 

the  left  hand  side  of  Eq  (3.16)  is  shown  to  equal 

r  d /'  [6  (f  +  /')  +  «(/-  /')]  (f\F)  =  Wtlfi  (/,  F)  .  (3.18) 

Jo 


Finally,  the  relation  between  W^(r,  /)  and  F ^  (kt )  can  be  expressed  as: 


/•OO  /*oc 

Wr4*(/.r)  =  4x  /  d ktktF[fi(k()  /  dr  cos(27r/r)  x 

Jo  Jo 

[2J0  (Artur)  -  Ja  ( kt  |f -  fir|)  -  J0  ( k ,  |r  +  ur|)] 


(3.19) 
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where  the  evenness  of  the  argument  in  r  has  been  exploited  to  simplify  the  integral.  Although  the  derivation  of 
the  phase  difference  spectrum  is  not  provided  in  their  paper.  GF’s  approximation  of  the  above  spectrum  can  be 
realized  by  setting  9  to  zero  and  performing  the  integrations.  Setting  9  =  0  is  equivalent  to  assuming  that  r  and 
v  are  colinear  over  the  integration  of  Eq  (2.1 1).  It  should  be  clear  at  this  point  that  such  an  assumption  is  clearly 
invalid. 

Having  determined  the  proper  form  of  1V^(/).  we  can  proceed  to  evaluate  the  single-sided  piston  and 
tilt  removed  spectrum,  W$(f),  using  the  GF  method  outlined  in  Chapter  II.  Recall  Eq  (2.11): 

W$U)  =  ~\j dr  Wtv  (f,  /)  (r) .  (3.20) 


Substituting  in  given  in  Eq  (3.19),  the  complete  integral  form  of  the  piston  and  tilt  removed  spectrum  is 


(/)  =  J  dr  4n  dkt 


.033C*kt2icPL 


,°° 

/  dr  cos  (2 x/r)  x 
Jo 


J  J  Jo  (Icf  -f  k~)  6  Jo 
[2 J0  (hvr)  -  J0  (kt  |r  —  vt\)  -  J0  (kt  |f +  rr|)]  (f) . 


By  making  some  convenient  variable  substitutions  and  interchanging  the  order  of  integration  we  arrive  at 
a  form  which  simplifies  the  numerical  integration  process: 


WiU)  =  3.040=  (f)(g)-!x 
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where /0  =  L2-.k  =  Dkt,t  =  =  jf,  and  C2  =  By  identifying  the  Fourier  integral  in  Eq  (3.22), 


we  can  also  define  the  oiston  and  tilt  removed  covariance  as: 


r°°  ke  {i&y 

d k - jT-  x 

.2\  6 


Cf(r)  =  1.54D2(^)"f  f'dxx  [  A6  Tv{x,6)  ! 

x  u '  Jo  Jo  Jo  /  /  „  \  i 

r+(^r) 

Jj0  (k\/  x2  +  t2  +  2 xt  cos  6  j  +  J0  {k\J  x2  +  t2  —  “Ixt  cos  0^  —  2  J0  (&<)j  (3.23) 


where  t  =  ^ .  Recall,  we  have  made  the  following  assumptions  to  arrive  at  these  equations:  plane  wave  input, 
frozen  flow  turbulence,  and  single  layer  atmospheric  turbu lance.  The  form  of  this  spectrum  and  its  numerical 
evaluation  will  be  discussed  in  Chapter  IV. 


3  3  Temporal  equivalent  spatial  filter 

Recall  from  Chapter  I  our  model  of  the  adaptive  optical  system  mean  square  error. 

0-2  =  f°°  d/  |1  —  (3.24) 

J  —  OO 

where \  V{  f)  is  the  temporal  power  spectrum  of  the  wave  front  phase  fluctuations  due  to  atmospheric  turbulence 
and  H  (/)  is  the  transer  function  of  the  adaptive  optical  system.  In  sections  3.1  and  3.2,  we  derived  the  temporal 
wave  front  spectra  W$(f)  and  W^(f).  Both  of  the  temporal  spectra  can  be  used  in  Eq  (3.24)  to  determine  the 
residual  wave  front  error.  We  can  determine  the  effect  of  aperture  piston  and  tilt  removal  on  cr2  by  comparing  the 
results  of  Eq  (3.24)  evaluated  for  W^(f)  and  W$(/).  In  order  to  make  this  calculation  we  also  need  to  define 
the  system  transfer  function  H(f  ). 

Notice  that  Wv(f).  W$(f),  and  //(/)  are  all  functions  of  the  temporal  frequency  variable  /.  We  define 
H(f)  =  Hs(f)  Hj(f ).  where,  as  stated  in  Chapter  L  Hs(f)  accounts  for  the  spatial  frequency  response 
of  the  deformable  mirror  and  Hr(f)  is  the  temporal  transfer  function.  We  have  assumed  that  the  spatial  and 


3-7 


temporal  transfer  functions  of  the  system  are  independent.  This  is  reasonable  if  the  computation  time  of  the 
actuator  weighting  functions  is  insignificant  compared  to  the  integration  time  of  the  wave  front  sensor. 


In  order  to  come  up  with  an  expression  for  Hs(f)  based  on  the  spatial  frequency  response  of  the 
deformable  minor  we  define  a  spatial  transfer  function  Hs{kt)  that  accounts  for  the  spatial  wave  front  correction 
characteristics  of  the  system.  In  section  3.1  we  derived  the  temporal  spectrum.  li>(/),  from  the  spatial  spectrum 
f\fl(kt)  using  relations  developed  in  Appendix  A  (see  Eq  (3.5) ).  Using  these  same  relations,  it  is  possible  to 
spatially  filter  ( k  < )  to  account  for  the  spatial  response  characteristics  of  the  defonnable  mirror  and  then  convert 
the  result  to  a  temporal  equivalent  spectrum,  W  s  (/)  ■  fit  other  words,  the  spatial  spectra  of  the  defonnable  mirror 
surface  is  modeled  as  Fs{kt)  =  F^(kt)  ■  Hs(kt).  Recall  that.  Eq  (3.24)  defines  the  spectral  power  difference 
between  the  incident  wave  front  phase  spectrum  and  the  spectrum  of  corrector  motion  of  the  deformable  mirror 
surface  (see  Fig  1.2).  Using  Eq  (3.5),  we  can  write  the  temporal  equivalent  spectrum,  Ws(f)  as: 

HW>  =  T  f  (fi+Wf)  -  <3'25) 

Finally,  the  temporal  equivalent  filter  Hs{f)  is  equated  to: 


Hs(f)  = 


(3.26) 


This  approach  is  relatively  simple  if  we  know  the  form  of  the  two  dimensional  spatial  phase  spectrum. 
F{kt).  Eq(3.4)  describes  Fv(kt),  the  spatial  frequency  equivalent  of  ( /)  Unfortunately,  we  do  not  yet  have 
a  spatial  frequency  representation.  F$(kt).  for  the  aperture  piston  and  tilt  removed  spectrum  W^(/).  Because 
W#  ( / )  varies  for  each  point  on  the  aperture,  the  random  process  described  by  ( / )  is  not  WSS  in  the  spatial 
domain,  and  as  such.  F^(kt )  does  not  exist.  We  recognize  the  requirement  of  stationarity  when  using  linear 
systems  approaches,  and  therefore  assume  local  stationarity  over  a  subaperture  centered  at  the  point  of  interest. 
This  assumption  approaches  validity  as  the  subaperture  size  decreases  and  is  necessary  if  a  linear  systems 
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approach  in  the  frequency  domain  is  to  be  used.  Most  important,  this  assumption  allows  us  to  hypothesize  the 
exitance  of  the  aperture  piston  and  tilt  removed  spatial  frequency  dependent  spectrum  ( k, ) . 


In  order  to  determine  the  effect  of  Hs(kt)  on  a  piston  and  tilt  removed  spectrum,  we  must  first  determine 
To  this  point  we  ha’  e  been  using  Eq  (3.5)  to  transform  a  spatially  dependent  spectrum  to  a  temporally 
dependent  spectrum.  Here,  we  postulate  that  it  is  possible  to  perform  the  inverse  operation  and  determine  F$(kt) 
given  W^(f).  In  Appendix  A,  we  determine  the  inverse  of  Eq  (3.5)  to  be: 

F<‘'»  =  -/"•••  a27> 

With  this  last  relation,  we  can  determine  F^(kt),  and  proceed  to  spatially  filter  this  spectrum  for  a  given  system 
spatial  response,  Fs  ( kt )  =  F^(kt)  ■  Hs(kt).  Eqs  (3,25)  and  (3.26)  can  then  be  used  to  determine  the  temporal 
equivalent  filter  Hs(f).  and  finally  the  residual  error  of  our  adaptive  optical  system  is  found  by  applying  Eq 
(3.24). 

3.4  Summary 

In  this  chapter  we  developed  the  mathematical  methods  needed  to  model  the  performance  of  an  adaptive 
optical  system  in  the  frequency  domain.  Section  3.1  developed  a  description  of  the  spectrum  of  atmospheric 
turbulence  induced  phase  fluctuations.  Section  3.2  developed  an  integral  representation  of  the  same  spectrum  that 
accounted  for  piston  and  tilt  removal  across  an  aperture.  Finally,  in  section  3.3,  we  developed  the  mathematical 
transformations  necessary  to  account  for  the  temporal  equivalent  spatial  response  of  a  system.  Fig  3.1  descibes 
the  process  developed  in  this  chapter  pictorally.  In  the  following  chapter,  these  methods  are  used  to  determine 
system  performance  for  several  spatial  and  temporal  transfer  functions. 
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Flow  diagram  of  the  frequency  domain  approach  used  to  determine  system  performance,  (a)  2-D 
spatial  spectrum  to  equivalent  temporal  spectrum,  Eq  (3.5);  (b)  2  D  spatial  spectrum  to  phase 
diffrence  spectrum.  Eq  (3.19);  (c)  GF  transform  accounting  for  piston  and  tilt  removal,  Eq  (2.1 1); 
(d)  determination  of  the  locally  stationary  piston  and  tilt  removed  spatial  spectrum,  Eq  (3.27);  (e) 
Fs(kt)  =  Fp(kt)  ■  (0  2-D  filtered  spatial  spectrum  to  equivalent  temporal  spectrum. 

Eq  (3.5);  (g)  determination  of  the  temporal  equivalent  spatial  transfer  function,  Eq  (2.11);  (h) 
spectrum  of  incident  wavefront  phase  (i)  determination  of  residual  error,  Eq  (3.24). 
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IV.  Numerical  Results 


In  this  chapter  we  apply  the  methods  developed  in  the  previous  chapters  to  determine  the  performance 
of  a  theoretical  adaptive  optical  system.  Before  Staining  results  for  system  performance  however,  it  is  still 
necessary  to  determine  numerically  the  form  of  W$(f).  Recall  that  W^(f  )  is  the  piston  and  tilt  removed 
temporal  spectrum  for  a  single  point  in  the  aperture.  These  calculations  are  very  time  intensive  and  only  few 
representative  points  on  the  aperture  will  be  evaluated.  Section  4. 1  details  these  calculations  and  presents  plots  of 
W${f).  In  section  4.2,  we  assime  forms  for  both  Hs{kt)  and  Recall.  Hs(kt)  accounts  for  the  spatial 

response  characteristics  of  the  deformable  mirror  and  Hj  ( /)  for  the  temporal  characteristics  of  the  control  loop. 
In  section  4.3.  we  evaluate  Eq  (3.24)  to  determine  the  residual  wave  front  error  as  a  function  of  system  response 
time.  These  results  are  presented  along  side  those  of  Fried  and  Karr  for  comparison  (2)(1 1). 

4  1  Calculation  of  the  piston  and  tilt  removed  spectrum 

To  numerically  evaluate  the  piston  and  tilt  removed  covariances  and  spectra.  C^(r)  and  for 

single  points  across  an  aperture,  we  must  define  the  following  constants:  C\ ,  C2 .  and  r  1 .  Recall,  Cj  =  .  and 

C2  =  7^  and  r  1  describes  the  location  of  the  subaperture  in  the  overall  aperture.  Because  of  the  difficulty  and 
time  involved  in  evaluating  Eqs  (3.22)  and  (3.23),  we  define  a  limited  number  of  rj  positions  on  the  aperture. 
We  will  evaluate  only  the  covariances  and  spectra  for  ^  =  0. 0,0. 225, andO. 45  where  the  radial  distanceshave 
been  chosen  to  lie  in  a  line  along  the  same  vector  direction  as  the  wind.  It  is  expected  that,  the  power  of  the 
temporal  wave  front  fluctuations  will  be  at  a  maximum  along  this  line. 

Although  many  different  values  are  quoted  for  the  size  of  the  outer  scale,  Coulman  et.  al.  find  that  L0 
has  an  upper  limit  of  about  5  meters  for  problems  associated  with  astronomical  seeing  (1).  From  Goodman,  the 
inner  scale  size,  l0 ,  is  on  the  order  of  about  a  millimeter  (5).  Using  this  information,  we  determine  C2  to  be  on 
the  order  of  4000.  The  spectral  power  at  frequencies  on  the  order  of  ( /„ )  “ 1  is  very  small  making  our  analysis 
insensitive  to  small  error  in  estimating  C2.  This  ratio  will  thus  be  held  constant  throughout  our  calculations. 
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Fig  4. 1  illustrates  the  effect  of  piston  and  tilt  removal  on  the  reference  input  spectrum.  Note  that  the  power 
in  the  piston  and  tilt  removed  spectra  increases  as  the  magnitude  of  the  radial  position  r  i  increases.  This  is  due  to 
the  higher  variance  of  the  tilt  correction  as  the  measurment  is  taken  closer  to  the  edge  of  the  aperture.  Also  note 
that  piston  and  tilt  correction  is  effective  up  to  temporal  frequencies  on  the  order  of  jr .  For  temporal  frequencies 
greater  than  ^  the  curves  have  the  familiar  f~8/3  slopeand  the  input  and  piston  and  tilt  removed  spectra  have 
the  same  power.  Equivalently,  we  can  determine  that  piston  and  tilt  correction  is  effective  for  spatial  frequencies 
up  to  jj  using  the  fact  that  /  -  (see  Appendix  A). 

The  parameter  C\  =  was  chosen  to  equal  2ir.  This  value  was  chosen  so  that  D  would  be  on  the  order 
of  two-thirds  meters  for  L0  =  4.  This  has  significance  because  we  can  compare  our  theoretical  low  frequency 
data  to  that  collected  experimentally  by  Greenwood  for  an  aperture  of  D  =  0.6  meters  (7).  Greenwood’s 
xperimental  data  plot  depicts  the  same  flattening  in  the  low  frequency  portion  of  the  spectrum  as  shown  in  Fig 
4.1  for  n  =  0.0.  This  flattening  occurs  at  a  value  of  about  10-2  waves2 Hz’-1  experimentally  as  compared 
with  our  theoretical  value  of  0 .0 1 3  waves2 Hz” 1 .  This  agreement  between  the  experimental  and  theoretical  data 
gives  us  confidence  in  our  spectrum  modeling  methodology.  Comparing  the  spectra  of  Greenwood  and  Fried 
described  in  Chapter  n  to  those  of  Fig  4.1,  we  find  that  accounting  for  the  vectoral  nature  of  Wtv{r,  f)  has  a 
significant  effect  for  low  frequencies  and  and  none  for  frequencies  following  the  f~6!3  asymtote.  Recall  from 
Chapter  n,  the  GF  piston  and  tilt  removed  spectra  follow  an  f4^3  asymtote  for  temporal  frequencies  less  than  ± . 

In  Fig  4.2,  we  allow  the  aperture  size  to  increase  to  that  of  the  outer  scale  so  that  C\  =  =  l.  This 

figure  shows  that  aperture  piston  and  tilt  correction  are  not  as  effective  at  removing  phase  distortions  when  the 
aperture  size  becomes  large  relative  to  the  outer  scale.  When  the  period  of  a  particular  spatial  frequency  is  less 
than  the  aperture  size  it  is  unaffected  by  piston  and  tilt  correction.  Thus  as  the  aperture  grows  larger,  the  power  in 
correspondingly  lower  spatial  frequencies  remains  unaffected  by  aperture  piston  and  tilt  correction.  In  Fig  4.1, 
piston  and  tilt  correction  removes  over  90%  of  the  of  the  spectral  power  present  in  the  reference  wave  front  at 
r  i  =  0.0;  whereas,  only  33%  of  the  spectral  power  is  removed  in  Fig  4.2.  When  comparing  the  spectra  of  Figs 
4.1  and  4.2,  it  is  important  to  remember  that  the  power  is  scaled  by  the  size  of  the  aperture.  Taking  D  into 
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account,  the  input  spectra  of  both  figures  have  the  same  value.  Realizing  this,  it  is  clear  that  a  smaller  percentage 
of  the  input  power  spectra  is  removed  in  as  D  approaches  L0. 


4.2  System  performance 

In  this  section,  we  specify  the  transfer  functions,  //(/),  used  to  represent  the  the  spatio-temporal  response 
characteristics  of  several  adaptive  optical  systems.  Recall  that  H(f)  =  Hs(f)  •  #t(/)  and  that  Hs(f)  is 
determined  from  Hs{kt )  using  methods  outlined  in  Chap  term.  The  spatial  transfer  function.  Hs{kt),is  defined 
in  subsection  4.2.1  and  the  temporal  transfer  function,  in  subsection  4.2.2. 


4.2.1  Spatial  filter.  A  deformable  mirror’s  spatial  transfer  function,  Hs(kt).  can  be  characterized  by 
the  size  of  its  subapertures  for  a  segmented  type  of  device  or  by  the  distance  between  actuators  and  the 
actuator  influence  function  for  continuous  surface  type  device  (13)(19)(10).  In  this  investigation,  we  idealize 
the  performance  of  the  adaptive  optical  system’s  deformable  mirror  by  approximating  its  spatial  frequency 
performance  to  be: 


Hs(kt)  = 


1,  l*«l<3 


(4.1) 


(  0, 

where  d  is  the  diameter  of  the  deformable  minor’s  subapertures.  This  estimate  of  deformable  mirror  performance 
is  consistant  with  Nyquist  sampling  rates  and  has  been  chosen  far  simplicity.  Additinally,  Tyson  finds  that  spatial 
filtering  models  essentially  the  same  as  Eq  (4.1)  are  adequate  for  characterizing  the  performance  of  adaptive 
optical  systems  (19).  This  filter  choice  also  helps  us  avoid  the  issue  of  a  nonstationary  p(t).  Recall,  it  is  the 
inherent  piston  and  tilt  removal  of  the  adaptive  optical  system  that  causes  p{t)  to  be  nan-stationary.  Because 
piston  and  tilt  removal  only  affects  spatial  frequencies  below  kt  =  we  can  consider  those  above  kt  =  2* 
to  be  stationary.  Because  our  chosen  Hs  {kt)  only  affects  frequencies  greater  than  k,  —  we  consider  our 
linear  systems  approach  to  be  valid. 


At  this  point,  we  are  ready  to  determine  the  temporal  equivalent  transfer  function  Hs{f).  Using  the 
methods  described  in  section  3.3,  we  could  determine  F^kt)  and  F$(kt)  for  the  spectra  plotted  in  Fig  4.1, 
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Figure  4.2. 


Temporal  (A)  covariances,  ,  and  (B)  spectra.  and  W$ ,  of  wave  front  phase  fluctuations  for 
the  case  of  =  1.  The  piston  and  tilt  removed  spectrum,  has  been  evaluated  for  f\  =  0.0. 


Figure  4.3. 

Temporal  equivalent  spatially  filtered  spectrum  Ws(f)  depicting  the  effect  of  Hs(k,)  = 
rect(^)  for  ^  =  2.  5.  10,  and  20.  The  circled  error  is  the  difference  between  Ws(f) 

and  our  approximation  of  the  spatial  filtering  effect,  rect 


multiply  by  Hs  (*<).  and  determine  Hs(f)  using  Eqs  (3.25)  and  (3.26).  Instead,  we  choose  to  simplify  this 
procedure  and  model  the  temporal  equivalent  spatial  filtering  by 


I 


I/I  <  JU 
0,  l/l>&- 


(4.2) 


where  A'  =  §■.  The  spatial  cutoff  point  of  Hs  (/)  is  determined  from  /  =  4^-,  (established  in  Appendix  A). 
Fig  4.3  shows  this  to  be  a  reasonable  approximation  for  a  representative  number  of  ^  values. 


To  make  a  gross  estimate  of  the  maximum  error  involved  in  approximating  Hs(f)  by  Eq  (4.2),  we  note 
from  Appendix  A: 

the  kF  (k) 


W(f) 


4?r  r°° 

v  Jmn 


1/2  ’ 


(4.3) 


From  Eq  (4.3),  we  can  deduce  that  the  spectral  power  at  any  temporal  frequency  /  is  the  weighted  sum  of 
the  spatial  frequencies  equal  to  or  above  k(  =  Because  the  spatial  filter  described  by  Eq  (4.1)  rejects 
frequencies  above  some  kc  0  =  J .  all  temporal  frequencies  will  suffer  a  loss  in  spectral  power  associated  with 
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die  spatial  frequency  band  from  kc  0  to  infinity.  Additionally,  since  the  weighting  of  these  spatial  frequencies 

is  greatest  for  fc  0  .  we  can  quantify  the  approximation  error  over  all  temporal  frequencies,  from  0  to 

£ 

fco.asW(fco)  x  fco  .  FromHg4.1  and  Eq  (3.9),  we  determine  that  W(/c0  )  =  (fco  )  3.  A  conservative 
estimate  of  the  maximum  error  involved  in  assuming  Hs(f)  has  the  form  of  Eq  (4.2)  is  (fco_ )~ 3 . 

4.2.2  Temporal  filter.  To  characterize  the  temporal  response  of  our  theoretical  adaptive  optical  system, 
we  assume  that  the  dominant  delay  in  the  actuator  control  loop  is  the  integration  time  of  the  wave  front  sensor. 
An  intuitive  linear  systems  approximation  of  this  integration  process  is  given  by  the  following  temporal  impulse 
response: 

hT(t)  =  rect  •  (4-4) 

Where  r  is  the  average  delay  between  wave  front  sensing  and  actuator  positioning  and  2r  is  the  length  of  the 
wave  front  sensor  integration.  In  the  Fourier  domain,  we  can  write  Eq  (4.4)  as 

HT(f)  =  sinc(27r/r)e'j2,r/T.  (4.5) 

The  phase  delay  term  in  Hj(f)  accounts  for  the  decorrelation  between  the  sensed  wave  front  and  the  corrected 
wave  front;  while,  the  sine  term  limits  the  temporal  bandwidth  of  the  wave  front  sensor  simulating  the  smoothing 
effect  of  the  integration  process. 

4.3  Calculation  of  residual  error 

In  this  section,  we  will  use  Eq  (3.24)  to  model  the  residual  error  of  adaptive  optical  systems  descibed  by 
the  spatial  and  temporal  transfer  functions  of  section  4.2.  The  temporal  frequency  spectra  of  Hg  4.1  will  be  used 
as  the  input  spectra,  W(/),  to  Eq  (3.24).  These  spectra  include  three  that  account  for  aperture  piston  and  tilt 
removal.  for  points  located  at  radial  distances  r  i  =  0.0,  0.225D,  and  0.45  D  along  the  direction  of  v 

and  one  that  assumes  no  piston  or  tilt  correction,  All  plots  of  system  performance  depict  wave  front 

,  .4 

error  measured  in  waves  versus  time  delay  of  the  system.  The  wave  front  error  has  been  normalized  by  (  —  )  * 
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and  the  time  delay,  r,  has  been  normalized  by  §■.  We  are  now  ready  to  evaluate  Eq  (3.24)  numerically  and 
determine  the  performance  of  our  system. 


4.3.1  The  effect  of  phase  delay.  In  Fig  4.4,  H  (/)  has  been  simplified  to  account  only  for  the  temporal 
delay  of  the  system,  Hr(f)  =  e~i2*ir  and  \H(f)\  =  1.  This  models  the  case  of  perfect  spatial  wave  front 
correction,  d  =  0  far  the  deformable  mirror  subaperture  size,  and  neglects  the  effects  of  temporal  bandwidth 


limitations  caused  by  the  integration  process  of  the  wave  front  sensor.  We  model  this  simple  case  so  that  the 
individual  effects  of  spatial  and  temporal  bandwidth  limitations  and  temporal  delay  can  be  distinguished  in  later 
plots.  The  curves  of  Fig  4.4  plot  wave  front  error  for  five  different  cases.  The  three  lower  curves  are  for  the 
piston  and  tilt  removed  input  spectra  The  upper  two  curves  assume  no  piston  or  tilt  correction.  The 

lower  of  these  last  two  is  for  the  input  spectrum  W^(f)  derived  in  section  3.1.  Recall,  H rv{f)  was  derived 
assuming  a  von  Karmen  index  of  refraction  fluctuation  power  spectrum.  The  uppermost  curve  is  a  plot  of  the  FK 
result,  (see  Eq  (1.1).  FK  assumed  a  Kolmogorov  index  of  refraction  fluctuation  power  spectrum.  The  FK  result 


can  be  written  in  closed  form  as 


(4.6) 
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where  Eq  (1.1)  has  been  rewritten  in  terms  of  r0  .  For  a  plane  wave  propagating  through  a  single  turbulence 
layer,  the  Greenwood  frequency  in  Eq  (1.1)  can  be  expressed  as  fa  =  0.4277 (6)(20). 

We  can  make  several  observations  concerning  Fig  4.4.  First,  the  error  described  by  the  curves  for 
ri  =  0.0. 0.225D.  and  0.45D  and  the  reference  wave  front  seems  to  approach  some  limiting  error  as  the  time 
delay  increases.  This  is  due  to  the  sensed  wave  front  becoming  uncoreelated  with  the  incoming  wave  front  as 
the  time  between  sensing  and  correction  increases.  The  small  increase  in  error  as  delay  time  gets  longer  is  due 
to  the  finite  power,  (see  Fig  4.1 ),  at  very  low  frequencies  in  both  and  W^(/).  The  upper  limit  on  the 

error  is  equal  to  twice  the  power  under  the  respective  spectrum.  Second,  the  error  described  by  the  FK  result  is 
higher  than  that  for  the  piston  and  tilt  included  reference  wave  front,  W^(f).  This  is  due  to  the  FK  assumption 
of  a  Kolmogorov  spectrum  which  has  power  approaching  infinity  at  low  temporal  frequencies  whereas  our 
assumption  of  a  von  Karmen  spectrum  limits  the  power  in  Wv(f)  to  finite  values  at  all  frequencies.  The  large 
increase  in  error  as  delay  time  gets  longer  for  the  FK  result  is  due  to  the  infinite  power  at  very  low  frequencies 
of  the  Kolmogorov  spectrum. 

4.3.2  The  effect  of  limited  temporal  bandwidth .  In  Kg  4.5  we  include  the  effect  of  the  integration  process 
to  the  phase  delay  effects  depicted  in  Kg  4.4.  As  stated  previously,  the  integration  process  of  the  wave  front  sensor 
affects  the  bandwidth  of //j(/).  Kg  4.5  still  assumes  the  case  of  infinite  spatial  bandwidth,  Hs(f)  =  1.  As  the 
integration  time  of  the  wave  front  sensor  becomes  longer,  the  temporal  bandwidth  of  #t(/)  becomes  smaller. 
This  decrease  in  the  temporal  bandwidth  of  Hj{f)  is  best  evidenced  for  long  time  delays  and  is  negligble  for 
very  small  r  due  to  the  small  amount  of  power  in  Wv(f)  at  high  frequencies.  As  r  becomes  large,  the  rms 
error  of  Kg  4.4  is  equal  to  approximately  \fl.  times  the  rms  error  of  Kg  4.5.  This  difference  occurs  because  the 
sensed  wave  front  and  the  corrected  wave  front  become  uncorrelated  as  the  time  between  sensing  and  correction 
increases.  This  loss  of  correlation  for  large  r  results  in  a2  =  a2en ,ed  +  (r2correcUd  where  a2ensed  «  v2corrected 
for  \Ht(J)\  =  1.  For  the  case  of  the  bandwidth  limited  transfer  function,  |//j-(/)|  =  sinc(27r fr),  a2stnsed 
approaches  zero  for  large  r.  Notice  that  Figs  4.4  and  4.5  are  plotted  for  a  and  not  a2.  The  difference  between 
the  two  plots  for  a  specific  curve  will  thus  be  \/2. 
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Figure  4.5.  Residual  wave  front  error  for  system  temporal  response  Hj(f)  where \Hs(f)\  =  1. 

4.3.3  The  effect  of  su'  perturesiie.  Figs  4.6  and  4.7  account  for  the  spatial  bandwidth  limitations  of  the 
deformable  minor  in  addition  to  the  effects  of  the  system  temporal  response  described  in  subsections  4.3.1  and 
4.3.2.  The  temporal  equivalent  spatial  bandwidth  of  Hs(f)  is  related  to  subaperture  size  by  Eq  (4.2).  Rg  4.6 
depicts  Jie  effect  of  an  aperture  to  subaperture  size  ratio  of  ^  =  10  on  system  performance.  Comparing  the 
curves  of  Rg  4.6  with  those  of  Rg  4.5.  we  note  that  the  subaperture  size  of  the  deformable  mirror  has  little 
effect  on  the  system  for  large  r.  As  the  temporal  response  of  the  system  becomes  faster  though,  the  finite  size  of 
the  subapertures  becomes  the  limiting  factor  of  the  system’s  wave  front  correction  performance.  For  7  =  10. 
Rg  4.6  shows  the  scaled  error  in  waves  to  be  just  less  than  0.01.  Rg  4.7  shows  that  as  the  subaperture  size 
grows  larger  the  error  also  increases.  For  very  limited  correction,  §  =  2,  the  scaled  error  plotted  in  Rg  4.7  is 
approximately  0.033  waves. 

From  Rg4.7,  we  can  determine  the  slowest  system  response  time  that  takes  full  advantage  of  our  system's 
spatial  response.  This  point  will  be  located  at  the  point  on  the  curve  where  significant  reduction  in  residual  error 
is  no  longer  achieved  for  faster  response  time.  We  will  refer  to  this  point  as  Tde>ign.  Comparing  our  results 
with  those  of  Fried  and  Karr,  we  find  that  the  rdt,,gn  predicted  by  our  study  is  similiar  in  magnitude  to  that  of 
the  FK  result.  The  coincidence  of  our  result  and  that  of  Fried  and  Karr  becomes  closer  as  we  approach  the  re 
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Figure  4.6.  Residual  wave  front  error  for  system  response  H  (/)  for  %  =  10. 


assumption  of  7  =  °°-  The  largest  discrepancies  between  the  FK  result  and  the  results  of  this  investigation 
will  be  for  partial  correction  systems  where  ^  is  less  than  around  five.  Fig  4.7  shows  that  the  FK  predicted 
temporal  response  requirement  is  almost  double  our  predicted  requirement  of  Tdesign  «  for  ^  =  2.  Thus, 
for  partial  correction  systems,  the  FK  result  does  not  predict  temporal  response  requirements  as  effectively  as 
for  higher  order  correcting  systems. 

From  Figs  4.6  and  4.7.  we  also  determine  that  accounting  for  piston  and  tilt  removal  is  not  necessary 
if  we  are  interested  in  the  area  of  the  the  curves  near  Tdeiign.  Kg  4.7  shows  that  rdt>ign  will  always  be  less 
than  ~  _  Corresponding  values  of  r  in  Kg  4.6  show  the  difference  between  residual  error  for  Wv  and  for 

the  three  curves  representing  W$  is  very  small.  This  difference  is  greatest  for  the  W$  curve  representing  the 
center  of  the  aperture  but  would  be  less  if  we  looked  at  an  average  W$  over  the  entire  aperture.  Recall,  we  only 
evaluated  at  three  points  in  the  aperture,  and  that  although  the  point  at  (r}  =  0.45D,  0  =  0)  has  the  worst 
case  spectrum,  the  point  at  the  center,  n  =  0,  is  best  case. 
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Figure  4.7.  Residual  wave  front  error  due  to  system  response  at  n  =  0. 

4.4  Summary 

In  section  4.1  of  this  chapter  we  evaluated  Eq  (3.22)  and  determined  the  piston  and  tilt  removed  spectra 
for  several  values  of  the  system  charateristic  constants  Ci ,  Ci ,  and  ?i .  The  resulting  spectra  were  found  to  be  in 
agreement  with  experimentally  gathered  data,  and  also  reflected  expected  performance  trends  of  piston  and  tilt 
correction.  In  section  4.2  we  defined  the  spatial  and  temporal  response  characteristics  of  several  system  transfer 
functions.  Section  4.3  presented  the  rms  residual  error  estimates  for  the  previous  section’s  system  transfer 
functions.  These  estimates  of  system  performance  were  found  to  predict  a  maximum  control  loop  delay  not 
greater  than  twice  that  predicted  by  the  FK  result.  This  agreement  of  our  error  predictions  with  those  of  Fried 
and  Karr  as  well  as  the  reflection  of  expected  performance  trends  validates  our  results. 


-iv/D 
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V.  Conclusions 


The  primary  objective  of  this  investigation  was  to  develop  a  frequency  domain  approach,  using  linear 
systems  methods,  that  describes  performance  effects  of  an  adaptive  optical  system’s  temporal  response  taking 
into  account  aperture  piston  and  tilt  removal  and  spatial  bandwidth  limitations  due  to  finite  subaperture  size.  We 
modeled  system  performance  as  the  mean  square  residual  phase  error  of  the  corrected  wave  front  <p  c(x,  t ): 

/oo 

df  (5.1) 

-oo 

In  Chapter  HI,  we  developed  spectral  representations  for  the  wave  front  incident  on  the  deformable  mirror,  U’(/), 
and  also  for  the  system  transfer  function  H(f).  The  unique  aspect  of  this  model  is  the  relative  ease  with  which 
performance  characteristics  of  different  spatial  and  temporal  system  response  functions  can  be  investigated. 
Although  methods  were  developed  in  section  3.3  to  account  for  an  arbitrary  transfer  function  //(/),  section 
4.2. 1  introduced  idealized  temporal  equivalent  spatial  transfer  functions  that  greatly  simplified  the  numerical 
evaluation  process.  Using  these  spatial  filtering  concepts,  the  effects  of  limited  system  spatial  bandwidth 
can  determined  more  efficiently  than  with  those  of  GF.  The  simple  multiplication  involved  in  accounting  for 
the  temporal  characteristics  of  a  system  make  this  model  especially  useful  for  comparing  the  performance  of 
different  temporal  transfer  functions. 

The  secondary  objective  was  to  address  when  use  of  our  approach  yields  substantially  different  results 
than  those  of  Fried  and  Karr.  Section  4.3  shows  that  for  r  greater  than  ^  the  FK  result  does  not  provide 
adequate  modeling  of  system  performance.  However,  response  times  of  systems  that  take  full  advantage  of  of 
the  deformable  mirror  spatial  response  are  less  .  As  stated  in  section  4.3,  we  call  this  temporal  delay  time 
'''design-  For  a  system  with  minimum  spatial  response,  =  2,  the  difference  between  the  FK  result  and  that 
predicted  by  our  approach  for  rde„gn  is  less  than  a  factor  of  two.  This  agreement,  in  predicting  the  residual 
error  for  rdf,ign  of  our  approach  and  that  of  Fried  and  Karr  validates  our  results.  Although  the  FK  result  is  more 
easily  evaluated,  use  of  the  methods  developed  in  this  thesis  provide  more  insight  into  the  performance  trade-offs 
of  systems  having  different  temporal  and  spatial  characteristics. 
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Additionally,  we  found  that  our  approach  could  be  significantly  simplified  if  we  were  interested  in  response 
times  less  than  r  near  Tje,ign.  The  most  difficult  part  of  our  approach  was  developing  the  piston  and 
tilt  removed  spectra  and  a  method  to  spatially  filter  this  spatially  nonstationary  spectra,  and  in  section  4.3.3,  we 
found  that  piston  and  tilt  removal  had  little  effect  on  the  error  predictions  generated  by  our  approach  for  response 
times  less  than  Thus,  our  approach  becomes  numerically  more  appealing  while  still  providing  the  ability 
to  compare  different  temporal  and  spatial  parameters.  This  is  especially  relevant  for  the  case  of  large  apertures. 
As  found  in  section  4.1,  piston  and  tilt  removal  have  significantly  less  effect  on  the  spectra  of  phase  fluctuations 
foi  aperture  sizes  aproaching  the  outer  scale  size  of  the  atmosphere,  L  0 . 

In  summary,  the  objective  of  developing  a  frequency  domain  model  of  adaptive  optical  system  performance 
has  been  met,  and  although  the  FK  result  adequately  predicts  the  residual  error  for  time  delays  of  interest,  our 
approach  determines  the  value  of  the  maximum  time  delay  that  sill  takes  full  advantage  of  the  spatial  response 
of  the  deformable  mirror. 
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Appendix  A. 


This  appendix  develops  the  methodology  that  will  be  used  to  transform  a  given  temporal  spectrum  to  an 
equivalent  two  dimensional  spatial  frequency  spectrum  and  the  reverse.  Our  development  follows  Tatarski’s 
almost  identical  development  for  the  three  dimensional  spatial  spectrum  (18). 

We  begin  this  development  by  defining  the  two  dimensional  spatial  single  dimension  temporal  correlation 
function  by: 

B  (f,  r)  =  (p(x,t)<p(r  +  r,t  +  r))  (A.l) 

where^(f.t)  isa  WSS  spatio-temporal  random  field  and  f  =  (Ax,  Ay).  Using  Taylor’s  frozen  flow  assump¬ 
tion.  the  temporal  changes  in  the  field  (  c,  t )  can  be  represented  as  simple  spatial  shifts: 

ip(x,t  +  r)  =  <p(x  -  vr,t) .  (A.2) 

where  v  repesents  the  magnitude  and  direction  of  the  field  flow.  Substituting  Eq  (A.2)  into  Eq  (A.l): 

B  (f,r)  =  B  (r  -  vt)  (A.3) 

The  three  dimensional  spectrum  of  B  ( r  -  vr )  can  be  found  through  the  Fourier  transform  relation: 

u(/c, /)  =  ^  ^  IJJ  dr  dr  r+2lr^T*f?  (r  —  vr)  (A.4) 


Performing  the  spatial  transform. 


where  F  (k)  represents  the  two  dimensional  spatial  spectrum  of  the  field  <p(x,t)-  The  important  relation  between 
the  spatio-temporal  spectrum  and  the  spatial  spectrum  for  a  frozen  flow  field  is  found  to  be: 


u  (k,  f) 


F{k). 


By  writing  B(r,r)  in  terms  of  its  Fourier  transform, 

B  (r,  r)  =  J I J  (it  J), 

we  can  begin  to  develop  a  relationship  between  the  temporal  frequency  spectrum, 

W(/)  =  /  d  fe^TB(r), 

and  the  two  dimensional  spatial  spectrum,  F(k).  Letting  the  spatial  difference  vector  r  in  B  (r, 
the  temporal  covariance  function  can  be  written  as: 


B  (0,  r)  =  B  (t)  =  J  d /  e^'T  [JJdKu  («,  /)  , 

and  the  relation  between  the  temporal  spectrum  and  the  spatio-temporal  spectrum  is  found  to  be: 

w  (/)  =  J J  d ku(kJ). 

Substituting  Eq  (A.6)  into  Eq  (A.  10)  and  switching  to  a  polar  integral  representation: 

W(f)  =  J™  dKK  J\eb(f  +  ^pjF  (K). 
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(A.6) 

(A.7) 

(A.8) 

r)  go  to  zero. 

(A.9) 

(A.  10) 

(A.  11) 


Assuming  u  (x,  f)  =  u(k,  /),  and  utilizing  the  evenness  of  the  6  integration,  the  temporal  frequency  spectrum 
can  be  rewritten  as: 

W{f)  =  J^  dxxF(x)  2jf  d06(/  +  ^cos0)  .  (A.12) 


Making  a  change  in  variables. 


W(f) 


=  2  /  dxxF(x)  f  - 

Jo  J- 1(1  — 


dA  _  2 x  /  2tt/ 
(1  _  a2)1/2  kw  V  +  kv 


(A.  13) 


where  A  =  cos  6  and  the  relation  6  (f  +  *£■)  =  (a  +  )  has  been  used.  This  simplifies  to 


W 


</>-*£ 


Making  the  variable  substitution  x  —  k;  + 


W 


(/)=—  f°<kt 

«  7o 


dx  xF  (x) 


we  write  the  equivalent  form 


(A.  14) 


(A.15) 


Eqs  (A.14)  and  (A.15)  will  be  used  as  the  means  to  find  the  temporal  equivalent  spectrum  of  a  filtered  isotropic 
two  dimensional  spatial  spectrum  in  Chapter  HI. 

For  the  three  dimensional  spatial  spectrum,  Tartarski  finds  an  expression  equivalent  to  that  derived  in  Eqs 
(A.14)  and  (A.15)  for  a  two  dimensional  spatial  spectrum. 


(/)  =  —  r  dx«*W 

V  Jl’  l/l 


W 


(A.16) 


where  (x)  represents  an  isotropic  three  dimensional  spatial  spectrum.  Using  the  relation  (17) 


_d_ 

dw 


/•■M")  dF  del)-)  ,  ddn 

/  dxF(x,w)  =  /  dx  — — t-  F(<fo,w)  — 1  -  F(^i,w)  — — 


(A. 17) 
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Tatarski  finds  the  inverse  relation  to  be 


${«)  =  — 


8tt3/c 


W' 


(A.  18) 


Unfortunately,  Eq  (A.  17)  can  not  be  similarly  applied  to  Eq  (A.  14)  to  find  the  temporal  spectrum  equivalent 
of  the  isotropic  two  dimensional  spatial  spectrum.  This  is  because  of  the  more  complex  integral  argument  of  Eq 
(A.  14)  involving  the  temporal  frequency  variable  /.  Knowing  the  three  dimensional  spatial  spectrum  equivalent 
of  the  temporal  spectrum  from  Eq  (A.18),  a  relationship  between  $  (k)  and  F  (k)  is  developed  instead.  Given 
the  Fourier  transform  relation: 


r 


d k  cos  (/c  •  r)  $  (k) 


(A.  19) 


where  T(r)isa  three  dimensional  spatial  correlation  and  r  =  (Ax,  Ay,  As),  the  two  dimensional  correlation 
can  be  found  from: 


T  (Ax,  Ay.  0)  =  fl(Ax,  Ay)  =  JJ  d/cx  drcy  cos  (kzAx  +  nyAy)  JdKz${Z) 


(A.20) 


The  relationship  between  the  two  dimensional  and  three  dimensional  spatial  spectrum  is  found  to  be: 


F («x,  Kyi *) 


J  d/c,  $(«t,  Ky.Kj)- 


(A.21) 


Finally,  we  determine  the  inverse  relation  of  Eq  (A.15)  by  substituting  Eq  (A.18)  in  Eq  (A.21): 


(A.22) 


where  $  (k)  =  $  (\At  +  “>d  k?  =  k7x  +  k 

In  summary,  this  section  has  developed  the  mathematical  means  to  convert  a  temporal  spectrum  to  an 
equivalent  isotropic  two  dimensional  spatial  spectrum  through  the  application  of  Eqs  (A.18)  and  (A.21)  and  then 


revert  back  to  an  equivalent  temporal  spectrum  by  Eq  (A.  14).  The  necessary  assumption  for  these  transformations 
was  only  that  of  frozen  field  flow. 
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